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ABSTRACT 


Randomly oriented filamentary composites form an 
important group of materials, f They find wide applications in 

w"' y ' 

both common and advanced engineering applications. But 

\ V 

unfortunately only a limited investigations have been made in 
this area.' A computer simulation technique using finite 
elements has been described here for the micromechanistic study 
of this type of materials to describe their in^plane elastic 
properties. A flat rectangular model with rectangular cross- 
section under tensile loading on the two ends is studied. The 
fibers have been randomly oriented using a technique based on 
random number generation. 

<1 The technique makes it possible to get the detailed \V 

stress distribution in the random structure which ip necessary, / 

* \ ■ 

. . . N A \ 

for the study of failure mechanics .' vit is found from the microj- 

mechanistic study that the maximum stress concentration factor 

in the random structure is much higher than that for regularly 

spaced fibers under identical conditions. However, the maximum 

shear stress seems to be of the seme order. For the simulated 

random structure, an interesting in-plane bending effect is 

found which disappears slowly as the fibers tend to align in a 


particular direction. This needs thorough investigations, both 


theoretically and experimentally 

^The macroscopic responses of randomly oriented 


filamentary composites have also been studied for the influence 



of biasness, residual compressive stress and volume fraction. 
Residual compressive stress has almost no effect on the modului 
(Young's or Shear). Tho increase in volume fraction with 
increased cross-sectional area of the reinforcements has littlt 
effect on the modulus for random structure. It slightly 
increases with increased volume fraction. But the biasness 
has great influence on the modului of the composites. They 
rapidly increase as the fibers align more and more 
particular direction. ^ 


in a 
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NOMENCLATURE 

E Young's modulus in psi 

G Shear modulus in psi 

a Normal stress in x or y direction in psi 

x >y 

cr Shear stress in xy plane in psi 

xy 

{a} Stress vector namely c x , cr^ and 

u, v Displacement along x or y direction in inches 
{u} Displacement vector namely u and v 
q Nodal displacements' in inches 

{q} Nodal displacement vector 

{S} Surface force vector in psi 

{£} Strain vector, normal and shear strains 

dV Volume integral 

V Volume of the element 

dA Surface integral 

A,-, Portion of the elemental area, where displacement is 

prescribed 

n { Compliment. ary energy functional 

U Strain energy 

[k] Stiffness matrix of the element 

Ep Potential energy functional 

A-^ Portion of the elemental area where surface traction 

is prescribed 
{Q} Nodal load vector 
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I Prescribed surface "traction in psi 

E- Modulus of fiber reinforcement or the matrix 

f ,m 

V_P Fiber volume fraction 

v Poisson’s ratio 

Poisson’s ratio of fiber reinforcement or the matrix 

i ,m 

[K] Overall or Master Stiffness matrix of the structure 


Bar indicates average value. 



CHAPTER I 


INTRODUCTION 


1.1 Preamble: 

Making composites is an old art but new science. Early 
civilisations are to be explored to find the genesis of composite 
As in many other walks of life, even modern science wonders over 
the artistic treasures made with older composites. Typical 
examples may be cited as" 1 ' straw-reinforced brick (early Isralite 
civilisation), Incas' platinum dispersed gold and silver ornament 
(early American civilisation), Damascus (18^ Century) and Japans 
(lb^ 1 Century) swords, Indian flint lock gun (17^ Century). But 
arts alone cannot survive unless they are supplemented by science 
As the aerospace science and technology started acquiring momentu 
there was an acute shortage of material to meet the new needs, i 
a result of the search of new materials, many composites were 
discovered and the search is still going on. Since Sixties a 
conscious and deliberate applications of the basic philosophies 
of composites with the different materials have made the progress 
of this group of materials very fast. At the moment composites 
are extensively used on a commercial basis even for domestic 
purposes. They are being widely used in industries like aircraft 
and missiles, architecture and construction, transportation, 

machine tools, agriculture, consumer products and appliances — 

2 : 
to name a few. ‘ 
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1.2 The Perspective: 

The importance of composite material is well established. 
Here, it may be of interest to review it from the point of nations 
interest. It has been rightly pointed out in the Report on the 
Science and Technology Plan of HOST group on special materials'? 
that the most widely applied composite, namely the Fibrous or 
Filamentary Composite , has a number of attractive features in the 
perspective of our socioeconomic needs. Traditional craftsmanshi| 
existing in the country can be economically and intelligently used 
for meeting many of the country’s needs without going for much of 
sophistication. The report points out: ” FRP (Fiber Reinforced 
Plastics) technology has an employment potential at both skilled 
labour and scientific development levels. The primary market will 
be local but export possibilities must surely arise later. M 
Typical uses have been cited in agriculture, public health, housir 
and transport apart from other advanced applications in different 
industries. The major organisations involved in the development a 
applications of fibrous composites include V.S.S.C., N.A.L., 
C.G.C.R.I., B.A.R.C. etc. apart from the advanced educational 
institutes like I.I.T.s and I. I. Sc. In terms of money about 3*6 
crores of rupees (about one fourth of total budget, Rs 15.15 crores 
on special materials) has been estimated for the development and 
application of the materials and Rs 10 lacs for basic research for 
5 years and $ 3,26,000 from U.H.D.P. for development of newer 


fibers . 
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1.3 The Problem; 

Under the perspective it is decided to synthesise the 
national need and latest technology available in the country. 
Finally, the problem is chosen from the area of filamentary or 
fibrous composites which apart from advanced technological appli- 
cations in defence and aerospace industries, will find equal use ir 
more common industries and areas with added scientific confidence. 

A computer simulation technique has been developed to 
evaluate the in-plane elastic properties of filamentary composites 
where short, discontinuous fibers are randomly oriented. The 
analysis is restricted to plane stress problems within the elastic 
region. Finite Element Method is used to measure the structural 
response of the composite system. 

The three compelling reasons for the selection of the 
specific problem are: : 

(1) To supply the 'intrinsic' material properties of filamentary . 
composites as the starting point of a reliable optimal design ! 
for advanced space and military applications, 

(2) To supply manufacturer and/or designer a tool to select the 
starting materials for their structural components of more : 

i 

‘ :•$ 

common industries, and I 

I 

(3) To help the fracture mechanics man with detailed stress 
distributions within the structure for analysis for a safer 1 
design. 

A master program has been developed for analysing filamen-j 

I 

tary composites with upto 12 phases and their arbitrary distributii 
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The program can be used for continuous reinforcements also without 
any difficulty. In actual experiments, the possible combinations 
of the matrix and reinforcement are literally thousands and one sel 
of experimental results are not, in general, applicable for the 
second set. Considering the time necessary for making the 
experiment and the cost involved, it is rarely used as a starting 
point of any design. On the contrary, the computer simulation 
technique established here, serves as a very effective tool for 
designing the desired structural properties from a large number of 
possible combinations of materials. Also, since this method is 
essentially a micromechanics study, it has its inherent advantages 
over other analytical methods which are essentially macromechanics 
conclusions . 


« 
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CHAPTER II 
LITERATURE SURVEY 

2.1 Fiber Reinforced Composites: 

_ , 4 

Broutman defines the composite materials as the man-made 
group of materials obtained by three -dimensionally mixing two or 
more chemically distinct materials which gives the properties 
unattainable by the individual constituents acting alone and there 
is no significant chemical reaction among the constituents. The 
base material in which the other materials are dispersed is known 
as the matrix and the dispersed phase(s) is (are) known as 
reinforcement ( s ) . Depending upon the aspect ratio (length by 
diameter) of the reinforcement the composites are categorised. The 
fibrous or fiber reinforced or filamentary composites are those for 
which the aspect ratio is much larger than unity. Depending upon 
use, the fibers are made continuous or discontinuous in filamentary 
composites . 

To fabricate filamentary composites the strong fibers, 
continuous or discontinuous, aro randomly or preferentially aligned 
in the matrix while the latter is in molten state and the whole 
system is then cured under specific temperature and pressure. 
Special care is taken to drive away the air bubbles from the end 
product. Apart from different varieties of plastics, the metals 
are also occasionally used as the matrix material. The choice of 
fibers can be glass , boron or other advanced materials , metallic oi 
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inorganic or natural filers. Consequently the choice of filamentary 
or fibrous composites is almost -unlimited. 

In any filamentary composite the load is mostly taken up 
by the fibers. The matrix takes only a small part, if anything at 
all, of the total load. The main role of the matrix is to bind the 
fibers together and add stiffness and rigidity in transverse 
directions. Also it protects the fibers from environmental and 
frictional degradations. For discontinuous filamentary composites, 
it is the matrix through which the fibers pick up the load, they 
carry . 

The host of problems the designer of filamentary composites 
needs to solve include complex microstructure, inhomogeneous 
distribution of different phases, various shapes and sizes of the 
reinforcements, their alignment and continuity, little understood 
interface, possible existence of voids and presence of residual 
stresses introduced during curing. This battery of problems are 
good enough to force any designer to make a series of assumptions 
in solving even the ’simplest' problem in filamentary composite. 

Now, in addition, if the fibers are assumed to be randomly distri- 
buted in the plane of composites, it becomes almost the last straw 
on the 'designer's' back! However, these host of problems never 
mean that no satisfactory solutions for composite structures exist. 
It only means that the designer must have the foresight to 
appreciate the relative importance of the simplifications he is 
making . 
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The final test of any structural design is the experiment. 
But the choice for filamentary composites Being numerous, the 
experimental technique often become.* an uneconomic design procedure . 
Consequently, a number of investigators have come up with their 
approximate theories based on different assumptions. Though there 

c: r-7 

are a number of analyses available for regularly spaced fibers, “ 
the number is very few for randomly distributed filamentary 

Q 

composites. Cox investigated the effect of orientations of short 
discontinuous fibers on filamentary materials and showed that for 
planar matrix all possible types of elastic behaviour can be 
predicted by composition of four approximate sets of parallel 
fibers. Without any consideration of the matrix material, he 
concluded that tho isotropic properties for randomly oriented 
filamentary composites are given by , 


E = E 1]L /3 


and G = 3^/8 


( 1 ) 


where E-^ is Young's modulus of the composite when all the fibers 

g 

are parallel to loading direction and continuous. Bishop and 

Erenchel^ have also been reported to derive similar results for 

very low transverse and shear moduli.'. Tsai and Pagano 11 have founc 

the above isotropic properties through invariant properties of the 

12 

composites. Christensen and Waals derived the effective stiffness 
using an interesting model. However, they neglect the geometric 
details for macroscopic response. The common point of all the 
analyses is the assumption of homogeneous, anisotropic material 
instead of phasewise isotropic but non-homogene ous distribution of 
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materials. So far as macroscopic response of the system, namely 

the evaluation of elastic properties, are concerned the above 

assumption might not be a resi gross assumption 5 but it may 

definitely lead to a non-conservative estimate in the sense that th 

detailed stress distribution estimated may completely be different 

from the actual one. For this reason only micromechanics study 

where exact phase distribution is used, has become very popular in 

the field of composites. Unfortunately this micromechanics study 

13-15 

has been used only for regular array filamentary composites. 

In this analysis the composite is treated as non-homogene ous but 
phasewise isotropic medium with random fiber distribution. So the 
local variations in the stress pattern due to different orientation 
and/or spacing of the fibers can be confidently predicted and 
consequently the mechanisms of crack initiation and propagation can 
be better understood. It is on this understanding a micromechanics 
study is undertaken to put more confidence on tho values of in-plar 
elastic properties. 

It is quite understandable that the future of composite 
materials, like many others, does depend on the success of their us 
as commercial materials. Most of the commercial composites are mac 
of chopped strand mat. Thus in the end product we have randomly 
oriented fibers in the plane of the composites. The method 
developed here gives the fabricator a means to assure his customei 
about the performance of the product since he knows both the 
microscopic and macroscopic response of the system now. ; 
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Although a few results are available to evaluate the 
isotropic properties of randomly oriented filamentary composites, 
this micromechanics approach adds a new dimension to them. Tsai 
and Pagano"' - '*' have suggested to use these isotropic properties 
whenever 'designed anisotropy' cannot he exploited. They have 
concluded that this isotropic strength is the minimum one which can 
be achieved for any composite system. So this analysis can he 
equally used to estimate the minimum strengths of the composite witl 
any combination of the matrix and the fiber(s). 

The well-known Finite Element Methods are used for this 
micromechanics study. A brief description of the numerical 
technique is in the next Section. 

2.2 Finite Element Methods: 

Finite Element Methods are now extensively used in 
engineering practice . Quite a number of programs are available 
against payment. Also there are at least three programs"^ 
freely available for solving a large number of structural problems. 
In addition, in every advanced centre of education and research, 
individual researchers are making thoir own programs. The method 
is based upon idealising a continuum structure with infinite 
degrees of freedom by a finite number of small sub-regions or 
divisions interconnected by a definite number of joints or points 
(the so-called Nodal Points), each with a predefined number of 
degrees of freedom. These regions or divisions are known as 
Finite Elements. 
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lor detailed discussion of the finite element methods, the 

17 ifi 

reader is referred to standard text books. ’ However, for the 
sake of completeness , the derivations of rectangular hybrid element 
and constant stress bar element are given below. 

2.2.1 Pian’s Hybrid Pormulation 

In this method the stress distribution in the element is 
assumed in terms of m (=s 5 here) undetermined stress co-efficients 
{p} as follows 23-25 

{o} = [E] {p} (2) 

This assumed stress distribution must satisfy the equilibrium 
equations. Also compatible boundary displacements (u> are assumed. 

{u} = [L] {q} (3) 

♦ 

The surface force {S} , along the boundary of the element , can then 
be derived from (2) as 



{S} = [B] {P} 

(4) 

Prom the stress-strain 

relation 



b 

i — i 

« t 

fl 

r*-\ 

(x) 

(5) 

the following matrices 

can be introduced 


[H] = 

/ [p] t [N] [E] av 

V 

(6) 

[T] = 

/ [R] T [1] dA 

(7) 


The total complimentary energy functional can then be expressed as 

n c = | {f3} T [H] {0} - {|3} T [T] {q> (8) 

20 

The condition of minimum complimentary energy is assured if 
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3H C = o i.e. {g} = [H]" 1 [1] {!> (9) 

Substituting (9) in the expression of strain energy, it is found 
U = | {|3} T [H] m = | {« T [I] T fH] _1 [T] {« (10) 

This means that the element stiffness matrix is given by 

[k] = [T] T [H]" 1 [T] (ID 

The stress is calculated from (2), namely, 

{a} = [P] [H] -1 [T] {q> (12) 

In this analysis, for hybrid rectangular elements (fig. la) 
the equilibrating stress is assumed as: 

0 X = q + P 2 y i o y = + fi 4 x ana o xJ = (13) 

The compatibility of displacements along element boundary is 
maintained by choosing the displacements as 

u = a-, + a 0 x + a y + a.xy | 

5 \ (14) 

v = a 5 + a g x + a ? y + a Q xy J 

Following the above mathematical steps , both the stiffness and 
stress matrices can be calculated. For the sake of convenience, 
[P], [R], [N] , [ T] and [l] are reproduced in the Appendix A. 

Since the final stiffness and stress matrices involve the matrix 
operations of these five basic matrices, they are done by the 
computer and hence not reproduced here. 
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2.2.2 Constant Stress Bar Element 

In Displacement method, the ’exact’ displacement in the 
structure is approximated by various degree polynomials. Eor 
Constant Stress Triangles (CST) only upto linear terms are taken. 

In terms of known model displacements, the assumed displacements 

... 17,18 
are written as 3 

{u} = [C] {q} (15) 

Brom which the strain matrix can be calculated in the form 

{£} = [B] {q} (16) 

Using the stress-strain relation, it is found 

{o> = [E] [B] {q} (17) 

Assuming that no body forces act, the total potential energy is 
given by 

% = f J [3] T [E] IB] { q } aV - J (q} T [C] T {T> dA (18) 

V *1 

where I is the prescribed surface traction in the portion A^ of 
the element . 

20 

If the principle of minimum potential energy is now used, 
then it is found that 

[k] = / [B] T [E] [B] dT (19) 

V 

and {QJ = J [C] T {T> dA 

A 1 

The stress matrix can be calculated from (17). 


( 20 ) 



1 ? 

For constant stress bar elements, it is assumed that the 
bar takes the axial loads alone. The displacement along the 
length of the bar is assumed as 

u = + (q 2 - q 1 ) * J (21) 

From (18), the matrix [0] is calculated as 

[C] = i [-1 1] (22) 

For isotropic, homogeneous one -dimensional bar element [Ej equals 
to E. So the bar clement stiffness is calculated from (21) and 
(19) as 

AE 
I 

where A and L are the cross-sectional area and the length of the 
bar respectively (fig. lb). If the bar is inclined at an angle, 

IV 1 Q 

the standard transformation ’ matrix is to be used to calculate 
the deflection pattern of the structure. However, for the problem 
at hand, no transformation is needed. 

2.2.3 Assembly and Convergence 

To solve for the deflection of the overall structure, 
all the element stiffness matrices are to be assembled as per the 
nodes and corresponding nodal degrees of freedom. For the hybrid 
rectangular elements, the degrees of freedom per node are 2. 

But for bar elements, the allowed degree of freedom in the 
problem at hand is one. So to assemble a bar and a rectangular 
element along the common nodes, a number of zeros are inserted 


1 

-1 


-1 

1 


(23) 
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in the stiffness matrix of the bar element to account for the 
non-existing degree of freedom in the perpendicular direction of 
the bar. This is done automatically in the program. The scheme 
is shown in Appendix B with reference to fig. lc . 

It is noted that both types of elements are standard and 
well tested for convergence. s * c ’ Actually for that 
reason only, thej ere picked up. So no special treatment is 
needed for the purpose. It is also hoped that on assembling these 
two types of elements in a single program the convergence is 
not affected*, after all the total degrees of freedom per node is 
not changed. Due to the choice of the typical problem, the total 
number of elements in the overall structure becomes fixed. The 
number of elements cannot be reduced to keep the structure 
unchanged. On the other hand, the number of elements are too 
large to be increased further witnout crossing the available 
computer memory. So no direct proof of the assembled convergence 
for the specific problem could be done. The only comforting 
point is that the chosen elements are small enough and in addition, 
an indirect proof of convergence is given in Chapter V. Also , 
a number of standard problems with knoxm solutions have been 
tested with the developed program. The results give equally good, 
if not better, results for these problems. The details of the 
program and its performance are given in the next chapter. 
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CHAPTER III 

THE PROGRAM AJSD THE PROBLEM 

3.1 The Master Program: 

A general purpose program has been developed for structural 
analysis under static loading with a switching character of using 
either CST or hybrid or both types of elements. Eor solving the 

19 

large number of simultaneous algebraic equations, Wilson's BANSOL 
is used unchanged. Rectangular displacement elements are formed by 
static condensation from four equivalent constant stress triangles 
(CST) . Rectangular hybrid element s^ are used with first five terms 
using linear normal stress variation and constant shear stress. 
However, extension of the hybrid elements for higher order poly- 
nomials is possible without much trouble. 

The program is then used for solving the problem of pure 
bending of a beam giveh' in Ref. 18, pp. 165-167 (fig. 2). The 
results are obtained for both displacement and hybrid elements and 
compared with the exact solution in Tables 1 and 2. As it is 
obvious that the hybrid formulation gives quite exciting result for 
this particular problem. Almost exact values of displacements and 
stresses are reached even with such a limited number of elements. 

The time taken is also slightly more than that required by 
displacement method for same number of nodes. 

With tMse encouraging results, the basic problem of 
evaluating elastic properties by finite element method is taken up. 



16 


As another test check, Young's Modulus E and Poisson's ratii v 
we- 

±& evaluated by displacement and hybrid method for an isotropic, 
homogeneous material. Uniform tension is given at the ends of a 
typical test specimen of 2" length and a square cross-section of 
1" x 1". The response of the system is calculated by finite element 
methods. Actually it is the computer simulation of the tension test 
Changing the material property and aspect ratio (length by width), 
if necessary, the elastic properties can, thus, be calculated for 
any material, homogeneous or non-homogeneous , isotropic or ortho- 
tropic or anisotropic. The only requirement for evaluating materia] 
properties by this method is a general purpose program which has 
been successfully developed. Obviously, a good computer is needed 
for getting the results . 

In the computer simulation of the tension test, the uniform j 
tension is replaced by a number of concentrated loads at the nodal 
points, the value of which is calculated by consistent loadings. 

The thickness of individual materials and the specimen is always 
taken as unity. In the program, provision is made to use upto 
twelve different materials. The advantages of automatic node and 
element generations have also been used. The total number of 
nodes and/or elements does depend on computer memory. All the 
calculations of the present analyses have been done with the IBM 
computers of I.I.T., Kanpur (7044) and V.S.S.C., Trivandrum (360). 
The maximum number of nodes, each with two degrees of freedom, 
is kept as 400. In this program, arithmetic average has been used 
to calculate the average stress and average strain. After all, 



17 


Young's modulus is the manifestation of the macroscopic response 
of a system under static loading. 

For the isotropic, homogeneous material, 25 nodes and 16 
elements have been used to estimate the elastic properties. From 
the deformation pattern, it has become obvious that larger the 
nodal number, the better the average. So finer mesh always gives 
a better result until it converges with the exact one. This is 
particularly true for evaluating Poisson's ratio. The numerical 
results of the elastic properties are given in table 5» The inter- 
esting feature to note here is as follows. So far as elastic 
properties, in general, are concerned, both the displacement and 
hybrid methods give almost the same result. But, it is already 
pointed out that the detailed stress and deformation pattern are 
quite different for these two formulations. This implies that for 
elastic properties, the displacement method is also equally good. 
However, for micromechanics study it is always preferable to go for 
hybrid method. 

From table 5, it is apparent that the values of Young's 
modulus obtained from either formulation is very close with the 
exact one. However, Poisson's ratio varies slightly from the exact' 
value. A finer mesh, giving a better average, will obviously lead 
to the exact value. It may be noted that this property is rather 
difficult to evaluate even experimentally. Fortunately , the effect 
of Poisson's ratio is not very much sensitive for structural 
analysis'. 
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3.2 Continuously Reinforced Composites: 

The method of evaluating clastic properties "by this 
computer simulation technique is then extended for continuously 
reinforced composites. There are a large number of formulae 

available for the prediction of longitudinal and transverse Young's 

5 6 27-29 

moduli . 9 9 Obviously all these results have come up because 

exact values for these parameters cannot be predicted due to 
associated problems like complex microstructure, interface charac- 
teristics, voids, fiber shape, size, distribution, alignment and 
continuity. These difficulties have made the standardisation of 
experimental techniques also very difficult. Consequently, the 
experimental results show a certain amount of deviation as shown in 
fig. 3. So to measure/predict the elastic constants, both the 
theoreticians and experimentalists are forced to accept a number of 
assumptions. The assumptions made in this analysis include perfect 
bonding and absence of voids. The detailed assumptions and their 
implications are treated later. , 

The details of the fiber orientations and loading systems 
for longitudinally and transversely reinforced composites are shown 
in figs. 4a and 5a. A very crude model with nine nodes only is usee 
to evaluate the longitudinal and transverse elastic moduli of 
continuously reinforced composites. Even with this crude model, 
the agreement with other results is quite reasonable. The finite 
element mesh and deformation patterns for these two specimens are 
shown in figs. 4b and 5b. The results are summerised in table 4. ' 



19 


It is quite interesting to note again that Young's modulus obtained 
in two finite element methods are the same. 

It is generally accepted that Paul's prediction for 
longitudinal Young's modulus is perhaps the simplest and most widelj 
applied result in the field of composites. The slight^ higher value 
obtained in finite element method is due to the crude mesh. It is 
believed that this will converge to Paul's result if finer mesh is 
chosen. To this context, the result of Ref. 28 is rather doubtful. 
On the other hand, the transverse Young's modulus obtained by the 
finite element methods is very much acceptable. Considering the 

variations in the experimental results,"* the result can very well 

leemel j 

be demand&d as the correct one. 

With all these data generated and compared, the program and 

its applicability in the field of composites have been fully assured 

So the actual problem outlined in the introduction is being stated 

and systematically analysed now. ' 

3.3 The Problem and Solution Technique: j 

I 

In chopped strand mat short, discontinuous reinforcements j 

i 

are randomly distributed in the plane of the composite (fig. 5a). 1 

Also composites with natural fibers or low-cost-high— sale s-volume , 

| 

commercial composites, in general, often have randomly oriented 
in-plane fibers. The much publicised low-density high performance 
filamentary composites mostly owe their desired strength in the 
exactness of the knowledge of load distribution. On the other hand, 
if the loading history is not fully known, as is the case in many 
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real life applications , the anisotropic characteristic of these 

materials cannot be exploited. Actually in such cases, it is 

always preferable to start with isotropic properties of randomly 

oriented fibers for reliable, optimal design. This type of 

isotropicity has been recommended to introduce through the thickness 

for launch vehicle applications.^ Similar isotropicity can be 

achieved with four suitably oriented laminae (fig. 6b). Tsai and 
11 rbt. 

Pagano have suggested to use these isotropic properties as a 

Cj 

measure of minimum stiffness capability which can be confidently 
used whenever, the advantage of designed anisotropy cannot be 
exploited. 

To facilitate both the microscopic and macroscopic analyses 
of such an important class of composites , finite element method is 
chosen as the solution technique. Young's modulus and Poisson's 
ratio are calculated by simulating pure tension test of a randomly 
oriented filamentary composites in the computer and then taking the 
arithmetic average of the test specimen. Shear modulus is calculate 
using the standard formula for isotropic material 


G- = E/ 2(l + v ) 


( 24 ) 


'l r\eL I cetha <; 

The bar imgsfbie© the average value . 

To simulate a tensile test specimen of randomly oriented 
filamentary composites, the fibers are to be randomly distributed 
in the plane of a typical test specimen of matrix material. To 
orient the fibers randomly in the plane of the composite, random 
numbers are generated from the standard library program of the 
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computer. Equal probability for any orientation out of n possible 
orientations is chosen. So, the total span of the random numbers, 

O-VveL 

greater than zero ix> less than unity, is divided equally in groups, 
1, 2, .... n. As a particular number is generated, the correspon- 
ding group for the number is identified. Since every group 
corresponds to a particular orientation, the orientations of the 
fibers are systematically identified from the generated random 
numbers and accordingly placed in the matrix test specimen. In 
this particular analysis, only longitudinal and transverse 
orientations have been chosen. Also it is assumed in the analysis 
that the fibers act as bars. So in analysing the problem by finite 
element method, the fiber volume fraction is to be introduced 
through the cross-section and/or length of the short, discontinuous 
fibers . 

After careful consideration of the tensile test specimens 
used for conventional, plastic and composite materials , the 
dimension of the simulated test specimen is chosen as 8" x 2" x 1”. 
The thicknesses of fiber, matrix and test specimen are all unity. 

In addition, it is decided to keep the length of all the fibers 
equal and as unity. But to justify the assumption of bar element 
for the fibers, the fiber aspect ratio (length by width) should be 
much larger than unity. This implies that all the fibers are 
plates (fig. 7a). Bow, to solve the problem of fiber volume 
fraction, the test specimen is divided into sixteen blocks (fig. 7b) 
Each of these sixteen blocks is filled with fibers of identical 
orientation. The number of fibers in a block depends on the fiber 
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volume fraction and chosen fiber width. As for example, for a 
chosen width of 0.1'*, corresponding to fiber aspect ratio 10, the 
number of fibers in each block will be 3 for a fiber volume 
fraction 0.3. Thus a methodology is established for simulating 
randomly oriented filamentary composites with the help of a compute! 

Poliowing the methodology and using the simulated tension 
test both the microscopic and macroscopic responses of randomly 
oriented filamentary composites can be studied. It is also decided 
to study the effects of a few important parameters on the macro- 
scopic response of the material, measured through Young’s modulus. 
The parameters considered include randomness, biasness, volume 
fraction and residual stress. 

3.4 Assumptions; 

Identifying the problem, the following assumptions are made 
for the analysis; 

(1) Both the matrix and the fibers are isotropic, homogeneous 
materials in bulk and in-situ. 

(2) All the plate-like fibers are geometrically identical. 

(3) The composite is free from voids. 

(4) All the fibers are perfectly bonded initially and remain so 
during loading. 

(5) The composite is in a state of plane stress and obey Hooke's 
law. 

In addition, to solve the problem by finite element method, 
the following assumptions are also added. 
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(6) Effects of stress concentration, at the ends of the fibers, 
are not considered explicitly. 

(7) The fibers tako only axial loading and the axial load developed 
is constant all through the length. 

Since the composite is in a state of piano stress, the 
plate-like fibers do not disturb the in-plane stress distribution. 
So, oven with curious geometries, the fibers are representative. 

13 

Actually, analysis with similar fibers has been done successfully. 
The other assumptions regarding the composite have proved to be 
realistic many a times. 

Although the effect of stress concentration has not been 
considered explicitly, it is believed that the elements near the 
fiber tips have been chosen small enough to arrest the steep stress 
gradient effectively. The fibers donot have any stiffness in the 
transverse direction in bulk. So the fibers in-situ too can be 
treated as bar, particularly if the fiber aspect ratio is kept high. 
Eor discontinuous fibers, the axial stress is not actually constant 
in the fibers. It is constant beyond the critical length only. The 
constancy of axial stress implies an effectively longer fiber. 
However, this slight increment in effective fiber volume fraction 
is rather compensatory since the treatment of the fibers as lines 
in the finite element mesh effectively increases the volume of the 


matrix . 
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The results of this analysis are, now, being discussed 
in the next two chapters. Before studying the effects of 
different parameters on the performance of randomly oriented 
filamentary composites, the next chapter is completely devoted 
to the detailed study of a typical case. 



25 


CHAPTER IV 

MICE. OHE CHANI C S STUDY OF RANDOMLY ORIENTED 
FILAMENTARY COMPOSITES 


4.1 Introduction: 

If high, strength short discontinuous fibers are randomly 
distributed in a low strength molten matrix and the whole system 
is cured at a specific temperature and pressure , the resulting 
randomly oriented filamentary composite is isotropic material. 

This type of isotropic material is made, in particular, with natural 
fibers like coir, asbestos and sishal. In the commercial applicatr- 
ions of composites with chopped strand mat (CSM), the end products 
are similar isotropic materials. Thus this type of composites 
constitutes an important class of materials which need to be 
carefully analysed. Also it may be noted that whenever the 
designed anisotropy of the composites cannot be successfully 
utilised, the best solution is to use this type of isotropic 
composites. j 

The material properties of similar composites were first 
8 

evaluated by Cox. But no account was made for the matrix. So the' 
analysis applied to paper are not reliable for modern composites ! 
particularly when the matrix carries a part of the total load, 
Kfenchel^ is reported to have derived similar results for very 
low transverse and shear modului. The isotropic properties were 
derived by Tsai and Pagano 11 also using invariant properties. In 
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12 

a rather recent report Christensen and Waals have derived 
similar result based on a model of equivalent volume fraction for 
randomly oriented fibers. 

The one important common feature in all the above analyses 
is the derivation of elastic properties of the composite without 
considering the local inhomogeneity of randomly oriented short, 
discontinuous fibers. So far as the macroscopic response of the 
composite system is concerned, the above theories predict it within 
reasonable- limits. But to the authors best knowledge, no theory 
exists fee considers the local inhomogenieties . It is noted that the 
most widely used assumption of anisotropic'i^^homogcneity of the 
composites cannot predict tho local stress variations since the 
actual system is phasewise isotropic but distribution of the phases 
is non-homogene ous . But the local stress within a composite may 
vary widely depending upon the local inhomogeneity and it is the 
local stress variations which control the Initiation and propaga- 
tion of cracks causing the ultimate failure of the structure. 

4.2 Analysis and Assumptions: 

To accomplish the detailed study of the stress distribution, 

a micromechanics study is essential. Finite element method is used 

for the purpose. The matrix elements are simulated through five 

term Pian hybrid^* ^ rectangular elements. The fibers, assumed to 

17 "i l 

take only axial loads, are taken as constant stress bar elements. 

A general purpose program has been developed for solving structural 
problems with statical loadings by the use of hybrid or displacement 
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or "both, types of elements. 

For micromechanics study of two-phase randomly oriented 
short, discontinuous filamentary composites, a typical tension 
test specimen (aspect ratio = length/width =4), is subjected to 
pure tension and the corresponding response of the specimen is 
studied at different cross-sections. To simulate the random fibers 
with specific volume fraction, the test specimen is subdivided 
into 16 equal 1" x 1" cells (fig. 7). In each cell the fibers are 
unidirectional and they represent the exact volume fraction of the 
total system. To notch with the memory capacity and also to 
maintain the fiber volume fraction of 0.3, three fibers were chosen 
in each cell. For all the fibers, the geometry is identical. To 
simulate randomness, the random numbers are generated from the 
computer with uniform probability. Equally it can be extended for 
any other statistical distribution. For the sake of simplicity 
only two orientations of the fibers (longitudinal and transverse) 
are chosen. If memory permits,, four orientations (0°, 45°, 90° and 
135°) are sufficient to exactly simulate the randomness. However, 
the results seem to be equally good with two orientations only. 
Using again equal probability, the orientations of fibers in all 
the cells are docidod from the generated random numbers. The' final 
simulated structure of the 'isotropic system' is shown in xig. 8a. 
Since plane stress condition is assumed for the test specimen, tho 
geometry of the fiber, chosen as unit length, thin plates, does 
not really matter. The thickness of the test specimen is chosen 
as unity . 

Vs 

4 , 
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The isotropicity and homogeneity of both the matrix and 
fiber are assumed for bulk and in-situ conditions. Also the fibers 
are assumed to be perfectly bonded and there is no ■void in the 
composite. In addition, the fibers are supposed to take axial 
loads alone and the axial stress developed is constant all through 
the length. The effect of stress concentration at the end of the 
fibers is not treated explicitly. However, the finite elements 
near the fiber tips are chosen small enough to arrest the steep 
stress gradient effectively. The whole analysis is limited in the 
elastic region. 

The tensile test specimen is subjected to uniform tension. 
The uniformly distributed tensile stress is lumped at the nine modes 
of the final finite element mesh of the specimen by the use of 
consistent loading. The finite element mesh, the co-ordinate system 
and the consistent loading is shown in fig. 8b. 

4 . 3 Results ; 

The response of the test specimen under uniform tension 
of 4000 psi is calculated by the general purpose program developed. 
The material properties are given below 
Matrix (Epoxy Resin) 

C 

Young's Modulus = E^ = 0.4 x 10 psi and 

Poisson's Ratio = v m = 0.33 

Filaments (Glass Fiber) 

Young's Modulus = E^ = 11.8 x 10 psi and 

Poisson's Ratio = = 0.197 
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Cross-sectional area of each fiber = 0.1 sq. inch. 

It may be noted that the fiber depth being unity in each case , the 
fiber width is 0.1 inch for this case. 

Volume fractions of the filaments = V^ = 0.3 • 

4»3«1 Deflection Pattern 

The vertical displacement of the top and bottom nodes and 
the horizontal displacement of the side nodes are shown in fig. 9. 

As expected, the longitudinal deformation of the test specimen is 
not uniform. The deformation is less along all the nodes where 
fibers are placed longitudinally. In other nodes, matrix takes most 
of the load since the fibers are transversely placed. The matrix 
being a weaker material, the deflection is more along those nodes. 

The transverse deformation has a distinct bending pattern, 
eventhough the loading is pure tensile. To pinpoint the bending 
effect, the horizontal deformation of the side nodes is separated 
(fig. 9) from the deformed test specimen shown in fig. 8c. Obviousl; 
the coupling between pure stretching and bending is due to the local 
inhomogeneities. The inhomogeneity gives rise to the shift of the 
centroid from the plane of action of the resultant force. The shift 
of the centroid also varies from section to section as shown in 
fig. 10. Depending upon the shift of the centroid, an in-plane 
bending moment is produced which pushes the nodes horizontally. 

The direction depends on the sign of the bending moment which 
obviously depends on the nature of local inhomogeneity. Accordingly 
the deflection pattern may take any curious shape (fig. 8c) 
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depending on the random distribution of the filaments. Also it 
should be noted that the Poisson's effect plays an equal role to 
determine the final deformed structure of the 'test specimen'. 

The transverse fibers always pick up the transverse load, irrespec- 
tive of its source and thereby limits the absolute magnitude of 
horizontal deformation. It may be of interest to verify the 
in-plane bending effect for randomly oriented filamentary composites 
both qualitatively and quantitatively. 

4.3.2 Strains 

Pig. 8c shows the deformed structure of the test specimen 
simulated by the computer for a randomly oriented filamentary 
composite. The deformation is arithmatically averaged out (as shown 
by dotted line) to calculate the macroscopic response of the 
specimen. The deformation is magnified hundred times in the figure. 
For 8" long test specimen, the average longitudinal strain is found 
to be 0.426 x 10 for a uniform tension of 4000 psi. The detailed 
variation of longitudinal and transverse strain in the specimen is 
shown in fig. 11. To calculate longitudinal strain, different 
gauge lengths are chosen, each measuring from the top plane (y = 0) 
of the specimen. The longitudinal strain is pure stretching. The 
average longitudinal strain calculated from the graph, is expectedly 
almost identical with overall strain measured from average defor- 
mations of the whole test specimen. The variation is only about 
6 per cent. On the other hand, the transverse strain is influenced 
by in-plane bending. So simple arithmetic average might be 
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misleading to calculate the transverse strain. 

To take care of the bending effect in calculating the 
transverse strain, rather the Poisson's effect, ideally an integ- 
ration technique is to be used from point to point deflection. 
Instead the graphical method is used here. The transverse strain 
is calculated at different cross-sections corresponding to different 
values of -y and plotted. The average strain is then calculated 
graphically. At transverse sections, corresponding to -y = 1", 2", 
3" etc., the transverse strain is calculated here and plotted in 
fig. 11. From this graph average transverse strain is calculated. 

It may be noted that the two strains vary by an order of magnitude 
almost. From those strains, Poisson's ratio is calculated as 

v S ~ = Q- ( ^C ._ x , lcJ = 0.1462 * 0.15 (25) 

e 0.426 x 10 * 

«y 

There is a definite variation of either strain at different planes 
or gauge lengths, but nowhere the variation is alarmingly sharp. 

This assures that the variation is due to local inhomogeneity and 
there is no strain concentrations . 

4.3.3 Stresses 

In continuous filamentarjr composites, the stronger fibers 
always take most of the load. The role played by the weaker matrix 
is to bind the fibers together and protect them from environmental 
and frictional degradation of strength. The matrix also assures a 
transverse stiffness which the fibers alone cannot offer. For 
short discontinuous fibers , the normal stresses in the fibers are 
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developed through the shear stress of the matrix. Thus shear stress 
plays an important role for short discontinuous filamentary compo- 
sites. for randomly oriented filamentary composites, the fibers 
are effective, that is they carry most of the load, only if they aro 
placed along the direction of normal stress. If they are placed 
transversely with respect to the normal stress, they donot contribute 

Hsc- skarirwj of 

anything to sfcsasi- the normal stress. In such cases matrix need to 

carry the whole of the normal stress developed. 

fig. 12 shows a typical distribution of the normal stress, 

dy. , developed in the direction of loading. The stress developed is 

normalised with respect to the average stress applied. The cross- 

section is DD at a distance of 2.625". The different cross-sections 

mentioned herein are shown in fig. 9. It is clear from the figure 

that the stress in the matrix is hardly 10 per cent of the applied 

stress wherever the fibers are there to pick up the load. The stress 

in the fibers goes as high as 4 times the applied load. On the 

other hand, the fibers donot carry any load wherever the fibers arc 

transversely placed (right part of fig. 12). Consequently tho 

applied average stress is quickly developed in the matrix. The 

stress concentration effect is also quite clear from the figure. 

The maximum stress concentration factor is found to be 4-5. It is 

noted that this value is approximately double than that £ox of 

14 

regularly spaced fibers. 

The in-plane bending moment generated due to local 

inhomogeneity causes tension and compression along the length of 

the specimen. Consequently, a also varies along a longitudinal 

«y 
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plane corresponding to a fixed value of x, even if a constant 
tensile stress is applied, fig. 13 shows the variation of cr along 
two longitudinal planes. Thus it is noted that the variation Of 
the stress distribution is very much controlled by the bending 
stress, rather the nature of local inhomogeneity . 

The distribution of c is rather cirri ous . As mentioned 
already, the . distribution of a is controlled by Poisson's effect 
and the in-plane bending effect caused by local inhomogeneity. But 
the local inhomogeneity varies from plane to plane for a randomly 
oriented filamentary composite. Also the quantitative value of 
depends, in addition, on the overall distribution of other fibers. 
Pigs. 14-16 show the normalised stresses, and at different 

cross-sections. The stresses are normalised with respect to applied 
stresses and plotted as per cent values. As in case of a , o y is 
also mostly taken by the transversely oriented fibers. This means 

*i v> IBk 

that the transverse fibers, ineffective t=& take applied load, now 
takes a . Thus for randomly oriented fibers, the composite is 
never adversely weak in any direction. At least a particular 
percentage of fibers is always there to pick up the load coming 
from any direction 5 thus making the composite strong in every 
direction. Hence the randomly oriented filamentary composite can 
safely be used whenever the designed anisotropy cannot be exploited* 
On the other hand, however, since the matrix will have to carry a 
part of the load, the fullest exploitation of the fibers cannot be 
employed. So the ultimate choice of the composite depends on both 
'the magnitude and direction of the applied forces. 
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It is interesting to note the antisymmetry of a distri- 

button in the planes of D-D and G-G (fig. 16), BE and BP (fig. 15). 

The same antisymmetry is also present for the shear stresses in 

those planes. These antisymmetries reflect the opposite directions 

of the in-plane bending moments originated due to the shift of tho 

centroid of the inhomogeneous plane in opposite directions. This 

implies that from the nature of the distribution of c and a , 

x xy 

definite conclusions about the nature of the local inhomogeneity 
can be drawn. Otherwise, if the nature of the local inhomogeneity 
is fully known, the possible variation of c and o can be 

X XJ 

anticipated. Prom the analysis, it is found that maximum tensile 
stress in the transverse direction has developed in the matrix and 
equals to 40 per cent of the applied stress. On the other hand, 
the maximum compressive stress has developed in the fiber and 
equals to 145 per cent of the applied stress. It may be recalled 
that compressive stress is rather beneficial for better bonding of 
the fibers in the composite. 

It is assumed in the analysis that fibers take only axial 
load. But the matrix elements are allowed to develop the 
shear stress which is assumed to be constant. Consequently, an 
estimate of the maximum shear stress developed along the interface 
for the load transfer mechanism can be made which is valid only if 
the bond between the fiber and the matrix i"S perfect. The maximum 
shear stress developed in the composite is found to vary from -30 
to +40 per cent of the applied stress. Since the shear strength 
of any material is much lower, the shear stress developed in the 
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composite can cause concern. Thus, for randomly oriented filame- 
ntary composite, the shear strength may "become a controlling factor. 
Once a crack is developed due to shear failure, the crack may he 
propagated across the composite cross-section hy the transverse 
normal stress, a . This only implies that care inust he taken, hy 

Jh. 

complete stress analysis, in designing randomly oriented filamentary 
composite. A complete analysis alone can safeguard the applicability 
of this group of composites. The detailed stress analysis, thus, 
forms an important part in the design of this group of composites, 
in particular. 

4.3*4 Macroscopic Response 

The macroscopic response of the <fcomposite is finally 
evaluated hy arithmetical averaging process. The actual and 
averaged deformation pattern of the composite is shown in fig. 8c. 

The applied stress is 4000 psi. So the average Young’s modulus 
can he calculated. The average transverse strain is calculated 
from the graph. The calculated values of different parameters are 
given helow; 

Average longitudinal deformation = 0.03410893" (from analysis) 
Average longitudinal strain = 0.426362 x 10*"^ (calculated) 
Average applied stress = +4000 psi (applied) 

zr 

Average Young’s modulus = 0.94 x 10° psi (calculated) 

mmO 

Average transverse strain = 0.060 x 10 (from graph) 

Average Poisson’s ratio = 0.146 (calculated) 
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Average shear modulus is calculated as follows: 


Gr 


E 

2(1 + v) 


0.94 x 10 6 

2(1 + 0 . 146 ) 


0.41 x 10 6 


psi . 


Table 5 shows the predicted values for Glass-Epoxy system and 
compares with other available results. 

Unfortunately, no experimental results are available for 

1 p 

comparison. Christensen and Waals have compared their results • 
with experimental results (Fig. 17) for low fiber volume fractions. 
The material properties used are close to the one used in this 
particular analysis. But results are not available for comparison 
for 7^ = 0.3. However, from the trend in the figure shown, it 
seems that the difference between the predicted and the experi- 
mental result is quite large for higher 7 f . Also a comparison 
of the experimental results with the prediction of Tsai and 
Pagano is available as shown in fig. 18. The comparison is made 
for random short length boron filament reinforcement in epoxy 
matrix. For 7^ = 0.4 and low l/D ratio (= 10 here), the predicted 
values are definitely higher than the experimental one. Under 
the perspective, the results obtained by this method are believed 
to be equally, if not more, reliable than the existing ones. 
However, the real advantage of the method lies in detailed micro- 
mechanistic study. 
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4.4 Conclusions: 


be 

( 1 ) 


( 2 ) 


(5) 


f-re. sen.fr 

Prom the previous analysis, the following conclusions can 
summarised 

A methodology has been established to simulate the randomly 
oriented filamentary composites with the help of the computer 
and finite element method has been successfully used for the 
micromechanics study of a computer simulated test specimen. 

Pian hybrid method demands more extensive study for application 
in the field of composites. 

The simulated test specimen has shown a distinct effect of 
in-plane bending due to local inhomogeneity. The effect 
might be smoothened out if the number of blocks in the 
transverse plane could be increased substantially. However, 
it might be equally interesting to find the possible bending 
effect, its nature and magnitude, through careful experimental 
measurement. 

The effect of stress concentration is predominant for 

randomly oriented fibers. Por the glass epoxy system it 

goes as high as 4.5. Similarly the effect of shear is also 

very high for random distribution. The shear stress nay go 

as high as 40 per cent of the applied stress. So a detailed 

in 


microme 


chanics study is particularly helpful to use^ randomly 


oriented filamentary composites with confidence. 
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CHAPTER V 

PARAMETRIC STUDY OP RANDOMLY ORIENTED 
FILAMENTARY COMPOSITES 


5.1 Introduction: 

The future of composite materials, like many others, does 
not depend on some specific, scientific applications. Rather it 
does depend on the possibility of extensive, commercial applications 
of the material. Again for successful commercial application, the 
confidence about the material performance must be gained through 
systematic, scientific analysis. If the fibers are randomly 
oriented in a matrix, the end product will be isotropic and 
homogeneous. But, obviously, the strength of the composite will 
not be very high in a particular direction. On the contrary it 
will be equally strong in any direction. So, this isotropy is 
useful whenever the designed anisotropy cannot be used. Also, for 
a safer, optimised study these isotropic properties can be used as 
a starting point. In short, the randomly ■ oriented filamentary 
composites have an important role to play in both commercial 
utilisation and scientific investigation. 

But unf ortunately , very little has been done for this 

O *] A "1 *1 

group of materials. Cox, Krenchel, Tsai and Pagano, Christensen 
12 

and Waals have derived individually the different in-plane 
elastic properties of the composite. To assume phasewise isotropic 
but distributionwise non-homogene ous composite material as 



anisotropic, homogeneous material is the common feature in all 
the derivations. Although the methods lead to rather reliable 
results, the detailed stress distribution in the composite cannot 
be achieved by any of the methods . A microscopic approach has 
already been discussed for microscopic and macroscopic analyses of 
randomly oriented filamentary composites. In addition, to the 
author's best knowledge, no systematic analysis is available to 
assess the effects of a number of important parameters on the 
performance of this group of composites. Here an attempt is made, 
with the help micromechanics study, to study the effects of 
randomness, biasness, fiber volume fraction and residual stress on 
the in-plane elastic properties of the composite. Tsai and Pagano " 1 " 1 ’ , 

IP 

Christensen and Waals have studied the effects of volume fraction 
in their reports. 

5.2 Assumptions and Analysis: 

The tensile test specimens of randomly oriented filamentary 
composites is simulated by the computer and subjected to pure 
tension within the elastic range. The response of the specimen is 
obtained using finite element method. The matrix elements, capable ; 
of taking linearly varying normal stress and constant shear stress, 
are chosen from Plan's five term hybrid element. The fibers are ; 

I 

I 

assumed as constant stress bar elements which can take only axial o 

I 

load. Both the fiber and matrix are assumed to be isotropic, ; 

homogeneous material and are perfectly bonded. No void is assumed 
to exist. 
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It is already mentioned that the study is limited on the 
effects of randomness, biasness volume fraction and residual stress 
on the in- plane elastic properties of randomly oriented filamentary 
Composites. This group of composites are fabricated by randomly 
distributing short, discontinuous fibers in a layer of molten 
matrix and subsequently curing the whole system at a specific 
temperature and pressure. Since no control is made in distributing 
the fibers in the matrix, the fiber placement definitely varies 
from specimen to specimen. Similarly, in the simulated test 
specimen, the fiber orientation also may vary depending on the 
random number generated and/or there placement, longitudinally or 
transversely. But, true randomness always converges to an 
isotropic system. So, for a particular set of matrix and fiber, bh/ 
composite material should always be the same isotropic and 
homogeneous material provided the above mentioned assumptions are 
adhered to. 

In most of the cases the isotropic behaviour* of composites 
with random reinforcement is used whenever the loading direction is 
either unknown or arbitrary. On the contrary, if the loading 
direction is not arbitrary and also fully known, the designed 
anisotropy of the composites is fully exploited. It is of interest 
to investigate the intermediate phases between the isotropic and 
designed anisotropic phases. In terms of statistics, for the case 
of isotropic composites, the probability of fiber orientation is 

f 4 , 

equal in any direction. Bor designed anisotropy, the probability 
of particular fiber orientation( s ) is (are) specifically higher 
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than others. In this study, only two fiber orientations are used.* 
The random numbers generated, which lie between zero and unity, 
arc accordingly divided into equal or biased groups to represent 
the fiber orientations. For isotropic composites the probability of 
fibers in the direction of loading is used as 50 per cent. As the 
biasness increases, the probability increases. Here the probab- 
ilities used are 60, 70, 80 and 90 per cent. In addition, fibers 
oriented in the direction of or perpendicular to the direction of 
loading are also analysed to represent the extreme cases of 100 
and 0 per cent probability respectively. 

The effect of increasing fiber volume fraction can be 
simulated by increasing either the number or the cross-section of 
the fibers. But unfortunately, due to limitation of computer 
memory, the number of fibers could not be increased to increase fibe 
volume fraction. So the cross-sectional area of individual fibers 
has been increased and subsequently analysed to study the effect of 
increased volume fraction. It is noted that increase of volume 
fraction by this method is highly artificial, particularly when the 
fibers are assumed as bar elements capable of taking axial loads 
alone. However, it is hoped that this shows the trend of the effect 
of increased volume fraction. 

Due to differential thermal expansions of the matrix and 
the fiber, a residual thermal stress is always developed. It is 
of interest to study this effect on the performance of the 
composite. In this analysis, the fibers are assumed to act as bar 
elements which can take the axial stress alone. So the residual 
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stress is incorporated as the initial stress in the fibers. 

Tho same g3.ass-epoxy composite test specimen (8" x 2" x 1") 
is choson for the parametric study. The detailed response of the 
'test specimens' are now analysed systematically. 

5.3 Parametric Effects; 

5.3.1 Randomness 

True randomness always converges to a specific isotropic 
system of composites. The simulation technique should also equally 
converge to the sane composite for identical matrix and reinforcing 
materials provided all the assumptions are adhered to in the 
process of fabrication of the composite and in addition the random 
numbers implying random fibers are sufficiently large. The second 
condition is somewhat important in the contradicting contexts of 
convergency and computer memory. 

To identify random fiber orientations, a set of random 
numbers are generated, low filling the sixteen blocks of a test 
specimen of the matrix transversely and longitudinally with fiber 
of volume fraction 0.3, two randomly oriented filamentary composite 
test specimens aro made. The test specimens are then put under a 
■uniform tension and averaging the response of the test specimens, 
■the Young's modulus is calculated. 

Pig. 8a shows the test specimen simulated by the computer 
when the random fibers are placed transversely. Pigs. 8b and 8c 
shows the finite element idealisation and the response of the test 
specimen under uniform tension. The Young's modulus is calculated 
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as 9.4 x 10 for this specimen. The second test specimen and its 
response to a uniform tension of 2000 psi are shown in fig. 19 
when the random fibers aro placed longitudinally. The calculated 
Young’s modulus is 9,6 x 10 psi. In this particular case, all the 
elements have been taken as displacement rather CST elements. It 
may bo recalled that in either the displacement or the hybrid model, 
the calculated value of Young's modulus may differ slightly, if 
at all . 

Uoting the wide variations in the two test specimens, it nay 
safely be concluded that the sixteen blocks are sufficient for 
showing the convergence of the material properties, at least with 
this set of random numbers. It is noted in this context that no 
convergence of the material responso has been shown in the 
application of the finite element method. Actually no choice is 
loft for showing the formal convergence . The computer memory is 
the only reason for which finer mesh could not be taken. On the ; 

other hand, the structure of the test specimen is so complicated 
that coarser elements could also not bo taken. However, the 
comforting point has been that the elements are small enough. How, 
with this randomness, the convergence is proved indirectly. The 
wide variation in structures rules out the possibility that the 
convergence of the test specimens to identical isotropic material 
is coincident. Rather it proves that, as expected, the randomly 
oriented filamentary composites have identical isotropic properties 
irrespective of its origin provided the above assumptions are 
a.dhercd to. So, even if, the structural details may vary from 
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specimen to specimen, the end products are the same for random 
fibers . 

5.3.2 Biasness 

It is noted that although equal probability is assumed in 

aligning the fibers of the computer simulated test specimen, the 

number of fibers in the direction of loading is more than that in 

the transverse direction. It is also noted that the respective 

proportions are ^ and These figures are used in calculating 

the isotropic properties of randomly oriented filamentary composites 

_ 11 

as derived by Tsai and Pagano. Physically also it is quite 
understandable that effectively the load carrying capacity of 
randomly oriented fibers will be definitely more than — since all 
the fibers except the ones in transverse direction will carry a 
part of the load, the share being dependent on fiber orientation. 

How, instead of dividing the random numbers generated 
equally and then using then to find the orientation of the fibers 
in all the blooks (sixteen in total) sequentially, they can be 
divided into unequal sets and use similarly to find the new 
orientations of the fibers. The physical interpretation is that 
more fibers are being oriented in a particular direction, say the 
direction of loading. This method of biasing the fibers in the 
known direction of loading and thereby increasing the stiffness 
of the composite in that direction is termed as designed anisotropy. 
As the biasness is increased, the stiffness is also increased. It 
is of interest to see how the two are related. To the author's 
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best knowledge no such, relation exists. 

The basic idea is to find a possible correlation. So 
such a set of random number is chosen, which on equal division 
gives almost equal numbers, 9 and 7, corresponding to the two 
directions. Also in addition, deliberately more number of 
fibers are oriented in transverse direction rendering them 
ineffective. It is then termed as 50 per cent biasness. Since 
the random numbers generated are between zero and one, 50 per cent 
biasness means that the numbers less than 0.5 and numbers equal to 
or greater than 0.5 represent the two perpendicular directions of 
fiber orientations. In a similar way biasness of any other 
percentage can be defined. It is again noted that this artificial 
method is deliberately chosen just to simulate the biased composited 
and consequently try to find a correlation between biasness and 
stiffness of randomly oriented filamentary composites. 

Tig. 20a shows the simulated composite with 50 per cent 
biasness. Tig. 20b shows the deformation pattern of the composite 
test specimen when subjected to a uniform tension. Here all the 
nodes of loading are made free. Trom the deformation pattern, 
the Young's modulus is calculated. Tigs. 21 to 24 show the test 
specimens with varying biasness of 60, 70, 80 and 90 per cent. 

All the test specimens are subjected to a uniform tension. Here 
all the nodes of loading ere horizontally restricted. The 
corresponding deformation patterns are also shown. The 
deformations are magnified 100 times. The composite with all the 
fibers placed transversly to the loading direction and the 
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consequent deformation pattern under uniform tension is also 
shown in fig. This may be treated as 0 per cent biasness. It 

is noted that the fibor volume fraction is always maintained as 
50 per cent. The Young's modulus for all the biasness are 
calculated. Tor 100 per cent biasness (when all the fibers are 
in the direction of loading) , the modulus is computed from the 
rule of mixture. Computing the relative modulii for different 
biasness with respect to the modulus of 0 per cent biasness, a plot 
is made between the relative Young's modulus and biasness (fig. 

25). The interesting plot shows a regular trend of increasing 
modulus with increasing biasness. Also the rate of increment of 
the modulus increases as the biasness increases. Tho relative 
increase is almost nil for the first half of increase in biasness. 
On the other hand the increment is almost exponential for increasini 
biasness beyond 50 per cent. However, the curve is rather 
qualitative. 

Another interesting point should bo noted in the study of 
the effect of biasness. It is already mentioned that the local 
inhomogeneity gives rise to in-planc bending moment in the test 
specimen. But increase in biasness, in general, means a reduction 
In local inhomo gene it y. As the local inhomogoneity decreases the 
bending effect also reduces (figs. 20-24) and as such, there is 
no bending effect when the fibers are continuous and placed 
longitudinally or transversely (figs. 4 and 5). Thus increase in 
biasness decreases the in-plane bending effect in addition to the 
increase in Young's modulus in the direction of the biasness. 
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5.3*3 Volume Fraction 

It is well-known that volume fraction plays an important 
role in composites, in general. If the fibers are continuous 
and aligned in the direction of loading, the Young's modulus 
increases linearly with increasing volume fraction (rule of 
mixture). But the relation is not so simple for transverseljr 
placed filamentary composites. She effect of increasing volume 
fraction on Young's modulus has been calculated from the macro- 
scopic response of randomly oriented filamentary composites by 

11 1 ? 

Tsai and Pagano , Christensen and Waals . Their results have 

31 

been reproduced in figs. 17 and I81 It is noted that the slope 
of the curve of Young's modulus - vs - volume fraction is about 
half for the random fibers than that for continuous composites 
with fibers in the direction loading. 

The effect is investigated in this .analysis also. But it 
is already mentioned that for the problem of computer memory the 
increase in fiber volume fraction is artificially simulated by 
increasing the cross-sectional area of all the fibers. Obviously 
this artificial technique cannot be used for very high volume 
fractions. In those cases, the fibers should also be treated as 
rectangular element instead of bar elements . 

As the volume fraction is increased by increased cross*- 
•sectional area alone, the deformation pattern of the test 
specimens does not change . Only the average increase in length 
decrease with increase in volume fraction. As a re, suit, figs. 8a 



48 


and 8c are representative of the test specimens and their 
deformation pattern. 

Pig. 26 shows the effect of increasing volume fraction 
on Young's modulus of randomly oriented filamentary composites. 

The rate of increase in the modulus is rather low. It may he 
noted that corresponding increase for longitudinally continuous 
filamentary composite is 4.56 x 10 psi. This means that increase 
in filer volume fraction hy 0.4 (from 0.2 to 0.6), the increase 
in the modulus is about 11.4 times that of the matrix for 
glass-epoxy system when the fibers are continuous and they are 
in the direction of loading. But the increase in the modulus 
of randomly oriented filamentary composites for the same increase 
in fiber volume fraction is only 1.27 x 10 psi which is about 
0.3 times the matrix modulus. Obviously this does not tally 
with the prediction of Tsai and Pagano 1 "*" which predicts almost 
an order of magnitude higher increase . Detailed experimental 
verification alone can throw some light on this prediction. 

However, the artificial way of increasing the fiber volume 
fraction might be a possible reason for the discrepancy. 

5.3.4 Residual Stresses 

It is reported that the residual stresses donot have any 

l x 

severe effect on the uniaxial strength of cross-ply composites. '' 
But no literature is available about the Influence of residual 
stresses on randomly oriented filamentary composites. Particularly 
this might have significant influence from the viewpoint of 
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fracture mechanics. So, the study is carried out for this 
parametric effect too. The same test specimen, shown in fig. 

8a, is studied under uniform tension when the reinforcing bars 
are subjected to varying degrees of initial stress. The initial 
stresses are all chosen as compressive and four values are chosen, 
namely, -250 psi, -750 psi, -1000 psi and -2000 psi. Since the 
initial stresses are assumed to be introduced due to differential 
thermal expansions, they are introduced in all the bars irres- 
pective of their orientations. 

Interestingly, even with such a wide variation in 

residual stresses, the Young's modulus of the test specimen 

5 

remains almost the same. The value is about 9.4 x 10 psi in all 
the cases. However, there is a very slow but definite tendency of 
reduction in the maximum stress concentration factor with incroasii 
initial stress of compressive type . Table 6 shows the summary of 
the analyses. It is noted that the reduction in stress concent- 
ration is very small compared to the increase in initial stress. 

The maximum normalised stresses, cr and cr increase almost by 

■a. 

the same amount. On one hand, reduction in maximum stress 
concentration is helpful, even if the reduction is small. On the 
other hand, increase in maximum shear stress is detrimental, 

! 

particularly because shear plays a very important role for the 
composites, in general. This means that inadvertently introduced 
residual stresses does not help the fabricator in the ultimate 
analysis. Also deliberate introduction of residual stresses will 
not equally help, as is evidenced from the trend shown in table 6. 
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5.4 Conclusion: 

The methodology established earlier for the analysis of 
randomly oriented filamentary composites by computer simulation 
technique, adds the reliability and opens up the potentiality 
for parametric study. The convergence of the non -homogeneous 
structure to the isotropic system is really interesting, parti- 
cularly when the number of randomly oriented fibers are small. 
From the parametric study it appears that the effect of residual 
stress can be neglected for all practical purposes for this group 
of materials too. It adds no advantage to designer even in 
terms of fracture mechanics. It appears from the analysis that 
the attempt to increase the ^oung's modulus by using fibers of 
larger diameter has little success. The rate of increment of 
the modulus is very poor compared to the increase in fiber volume 
fraction. However, the modulus can be very readily increased 
if fibers are aligned in the direction of loading. A slight 
alignment improves the modulus to a considerable extent and the 
rate of increment follows a definite pattern. 

It may be noted in conclusion that the shear modulus also 
follows the same trend under the influence of the above 
parameters . 
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CHAPTER VI 

CONCLUSION'S 


6 . 1 Summing-up : 

Randomly oriented filamentary composites form an 
important group of composite materials. A methodology is 
presented here to make the micros true tural analysis and study 
the macroscopic response of these materials by computer simulation 
technique. It is found that the random fibers give rise to much 
higher stress concentration factor in the specimen than that if 
the fibers are regularly spaced. However, shear does not play 
any additional detrimental role in this case. The interesting 
feature revealed in the analysis is the in-plane bending effect. 

The local inhomogeneities in the simulated test specimen give 
rise to such bending effect. But in actual system the random 
fibers are thousands in number unlike two, chosen in the 
transverse direction of the idealised specimen. This may smoother, 
out the local inhomogeneities and the bending effect may disappear. 
This needs to be treated in further detail if a better computer 
with more memory than the one used is available. At the same time 
careful experimental investigation may also throw some light. 

The microscopic analysis presented here shall always give a 
better insight about the performance of the random structure. 

It is also noted that the same analysis technique— can— be eq ually 
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extended without any additional trouble for regularly spaced or 
continuously reinforced composites. 

In the macroscopic responses; under the influences of 
difierent parametric effects? it is noted that compressive 
residual stress has no effect and increased fiber volume fraction 
due to increased cross-sectional area lias little effect on the 
loung’s modulus of the composite. But biasness has got a very 
significant influence on the modulus. So little success in 
aligning the fibers in the direction of loading will enhance the 
modulus considerably. 

6.2 Recommendations: 

The success of the analysis can be still improved by 
incorporating the following simple modifications for future 
investigations . 

In the finite element analysis the hybrid method gives 
a definitely bettor estimate of the stress distribution than that 
given by the displacement? rather constant slfcress triangle? 
method, similar Ijr it can be seen if mixed method gives any 
better estimate. In this connection, it might be noted that in 
the mixed method both displacement compatibility and stress 
continuity is maintained between different phases of the 
composite which are supposed to play an important role for this 
group of materials. In hybrid method? however? stress continuity 
is sacrified to satisfy the equilibrium equations. In addition? 
instead of using constant stress bar elements as fibers? a 
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tetter bar element# can be used* particularly the one along 
which shear is allowed and/or variation of normal of stress is 
incorporated. 

Only with two fiber orientation, convergence is 
achieved in the analysis. But the recommended practice is to 
use four orientations. So it will be of interest to see the 
changes, if any, in the results when four fiber orientations 
like 0°, 45° > 90° and 135° are used. Obviously this will require 
triangular elements implying larger memory. 

Finally the methodology can be equally extended for 
evaluating the isotropic properties of randomly stacked laminated 
composite provided suitable changes are incorporated in the 


program 
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AP PENDIX - A 

The i oil owing aro the matrices used in deriving the 
stiffness matrix oi five term hybrid rectangular elements? 
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APPENDIX - 3 


Pi c . . 1c nb ows the assembly of a bar and a rectangular 
element, lor the bar IJ which can take only axial stresses, the 
stiffness matrix will be of the following form. It is noted that 
tile stiffness elements of the bar along 1 and 3 degrees of freedom 
are zero 
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But for the hybrid rectangular element IJKL, the stiffness elements 
exist for all the eight degrees of freedom. Consequently the 
stiffness matrix takes the form 
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In order to assemble , the elements of the stiffness matrix 
ox identical decrees of freedom are to he identified and then 
added to £et the overall stiffness matrix of the assembled 
structure. So for fig. lc, the overall stiffness matrix becomes: 
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Similarly a horizontal bar element like JK can also be 
assembled with the corresponding addition of stiffness elements 
along the common degrees of freedom of 3 and 5. In the program 
this assembly process is done by the computer. 
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•IG. i GENERALISED FORCES & DISPLACEMENTS OF cj. PIAN'S 
RECTANGULAR HYBRID ELEMENT b CONSTANT STRESS BAR 
iLEMENT & c THEIR ASSEMBLY. 
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FIGJ5q CONTINUOUSLY REINFORCED COMPOSITE 
WITH LONGITUDINAL FIBERS a Jb ITS DEFORMATION 
PATTERN. 
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FIG 7.0 TYPICAL FIBER GEOMETRY a _b_ THE MATRIX BLOCKS 
TO BE FILLED BY FIBERS. 





FIG- 8a -THE TEST SPECIMEN: COMPUTER 
SIMULATION OF RANDOMLY ORIENTED 
FILAMENTARY COMPOSITE WITH V f =307. 
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oriented filamentary composite test specimen. 

FGt.Sc. DEFORMATION PATTERN (AVERAGE & ACTUAL) 

OF THE TEST SPECIMEN (DEFORMATION MAGNIFICATIONS® 



HORIZONTAL. DEFLECTION 

PATTERN OF SIDE NODES D 
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HORIZONTAL DEFLECTION F 

PATTERN OF CENTRAL NODES G 



FIG. 9 SEPARATED BENDING DEFORf' 
BETWEEN PURE STRETCHING & BE? 
CROSS-SECTIONS USED LATER ON. 
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FIG. 11 VARIATION OF £„ AND £ y FOR THE -TEST SPECIMEN 
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FIG 14 DISTRIBUTION OF NORMALISED STRESSES AT D#E fiWf 
CRbSS'SECTIONS SHOWN IN FIG- 9 



- 4 - 











yl, 







\ 

j*C ,• 

* 

<T - ' 

1 

, 

s 

» \ 
'T 

A - • .’ I 

p 

! - . -.-*>4-. 

* A* * * ~ " 

. r* 

**• *«w> - . 

- -aV 

l 


/ V \ 


/ \ 



\ ' 






-~W— 




r; 


, f I? 

A 

v. 

» s \. 


#■ -♦ +— 

-V 

I 4,.' /’ - f'< J 


H-H 

-» IN. 


\ \\ 

* % 

\ \ s. 


7 

20 b 


\ 

1 \ ‘ 

'it / / 

\ \ d > 

\ N / ! 

\ \ > / 

\ > ' 

X / 

' ' \ / 

' / \ / 

\ / \ t 

\ / w 

0 ¥ 


Gy = +4 KPS! 
o D-D 

• G-G 

* H-H 

a I~i 

— <J X ‘ VARIATION 
CF XY VARIATION 
AVERAGE (J v 


FIG. 16 DISTRIBUTION OF NORMALISED STRESSES AT DIFFERENT 
CROSS-SECTION SHOWN IN FIG,' 9 




FIG. 17 COMPARISON OF EXPERIMENTAL 
RESULTS AND DIFFERENT PREDICTIONS 
FOR GLASS POLYSTYRENE SYSTEM. 
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FiG. 18 COMPARISON OF EXPERIMENTAL RESULTS a 
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INCREASING BIASNESS 

FIG. 25- EFFECT OF BIASNESS ON YOUNG'S 
MODULUS (V f =0.3) 
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FIBER VOLUME FRACTION 

FIG. 26 EFFECT OF FIBER VOLUME FRACTION 
ON YOUNG'S MODULUS FOR RANDOM 
FIBERS. (GLASS EPOXY SYSTEM) 



COMPARISON OF LISxI«*CBxIX r TS SELECTED NODE 



H 

j — j 


ft 

P 

pq 


Eh 


JO 

EH 

t~n 

H 

O 

ft 

P 

ft 

EH 

O 

ft 

ft 

ft 

CO 


CO 

r '~*l 

CQ 

CO 

ft 

n 

:o 


ft 

o 


ft 

o 

CO 


3 

(=~i 


o 

o 



p 




c 




♦H 



-P 



Cd 

H 



e 

$ 



p 

4 



o 



ft 

— 


I 

rH 



CD CP 
'H O 



CD <d 



P 

h_J0_^ 

— 

♦ 1 







ft 

1 


in 

M 

«> % 


nj 

b 

CQ 


• 


-P 


CM 


Pi 




0 


*o 


a 


m 


0 


o- 


H 


m 


p 


» 




in 


C H 


•v ✓ 

ft 

O 



b 

u 




0 


-p 


ft 


Pi 


g 


*r{ 

■ 



o 

ft 

, M 


*ri 


b 


0 



*<K 

ft 





ri 




*H 



Of 




m 

0 



0 

0 


ft 

ft 

0 


v W 

o 

0 


b 

ft 

0 

in 



H 

o 


ft 

-P 

* 


o 

CO 

o 


fH 


CN 


0 


in 


ft 


CM 


gj 


H 


ft 


H 




s. 

ft 




b 



p 




-p 




p{ 




*rl 




O 




ft 

K 




b 



0 




m 




cd 

1 



o 

i 







0 


0 






tS B 




ri 


•H H 




o 


ft 0 ft 



0 o 

EH 

rJ 2 +3 

O 



CQ ft 

uQ 

>5 

cd 



H -P 

O 

ft m 

X 



o cd 



ft 



PM P=! 





CO 







P 

0 

ft 

ft 

i 


O 

2 

rH 

iH 

5 


ft 

H 

0 

0 



ft 

2 

CQ 

CO 




ft 





2 

o 

in 




t *! 

PH 


1 



ft 



O 



ft 

0 


H 



a 

OQ 


’ft 

O 


A 

Pi 

CM 

M 0 



P 

p *H 

in 

m 



P 

O 0 

• 

in 





rH 

♦ 



ft 


rH 

1 

i — l 



ft 

l 



H 

1 o 





ft 

CQ *H 


m 



P 

ti -P 


i 



P 

n 3 o 


o 



H JP 


rH 



>H 

ft ft 


M 



P 

0 

in 

m 

o 



CD Ph 

00 

m 



P 

tiO 

H 

CO 



P 

a 5 0 


M 0 


H 

EH 

H 0 

# 

• 

, 

H 

A 

H 

<P 

CM 

CM 


H 

ft 

i> © 




B 

-4 t> 

M 0 

cn 

o 

m 

O 


o 

Ch 

o 

p 

ft 

i d 

• 

* 

• 


. 4 

*H *H 

H 

<n 

o 

O 

W cd 

G"> 

ft 

in 

EH 


Pi U 

MO 

o- 

O 


ca 

o -p 




p 

P ui 


in 



p 



i 



eh 

0 H 


o 



p 

two Cd 


i — ( 



p 

cd Pi 


{xj 



ft 

H *h 

m 

cn 



o 

0 ft 

M 0 

o 

O 


ft 

t> 2 


in 



Ph 

P -p 

o 

O" 




"’'“****' 

• 

' * 



o 

0 

CM 

ft 



M 

0 

i 



ft 

0 




CO 

Sh 


m 

i 



_ H 
■*-i 

-P 




p 

CO 


o 



p 

0 

£0*H 


rH 

W 



ft 

(T\ 

MO 



o 

cj 0 

in 


O 


ft 

Pi ft 

CM 

o 



0 

M 0 

H 



o 

> P} 

a 




CO 


m 

} 

ft 

1 



H 

P 

ft 


cn 

O 

O 


ft 

0 

ft 

o 

# 

o 

o 

* 


ft 

O 

H 

o 

o 


o 

o 

rH 

m 

in 

Lf\ 



P 

CM 

CM 

CM 



-P 

0 

a 



H 




H 

H 




H 

H 

H 





03 o 
m h- 

« • 

o o 


o 


in 

in 

o 


o 

H 

m 

rH 


O 


W 

rH 

X 

C- 

H 

CO 

b- 

H 

H 

O 

• 

• 

• 

m 

in 

in 


Cn 00 

mo cn 

CT> H 

in vo J 

O O I 

o o 

* ♦ 

o o 


CM 

m oo 

CM CM 

in in 

H- | 

in in i 
I — T H 

o o 

• « 

O O 


O 

O 

• 

o 

o 

o 

CO 


o 

o 


o 

o 

o 

CO 


1 

I 


-p 

rt 

CD 

& 

CD 

Q 

0$ 

H 

P4 

CQ 


V3 


u 

ft 

& 


-p 

o 


cd 

ft 


.35 




EFFECT OF COMPRESSIVE INITIAL STRESS ON THE YOUNG'S MODULUS 
OF RANDOMLY ORIENTED FILAMENT ORY COMPOSITES 
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